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Abstract. In this article, the continuous time random walk on the circle is studied. 
We derive the corresponding generalized master equation and discuss the effects of 
topology, especially important when Levy nights are allowed. Then, we work out 
the fluid limit equation, formulated in terms of the periodic version of the fractional 
Riemann-Liouvillc operators, for which we provide explicit expressions. Finally, we 
compute the propagator in some simple cases. The analysis presented herein should be 
relevant when investigating anomalous transport phenomena in systems with periodic 
dimensions. 
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1. Introduction 

Continuous Time Random Walks (CTRWs) ([1], [2]) are models describing the motion 
of individual particles (or any other conserved quantity) probabilistically. Since their 
introduction they have found many applications in physics, among other reasons, due to 
the fact that they allow exploring possibilities which go beyond the classical paradigm of 
diffusive transport. Their simplest realization, known as a separable CTRW, is defined 
in terms of two probability distribution functions (pdfs): the step-size pdf, p, and the 
waiting-time pdf, ip. It is well-known that when the step-size pdf is Gaussian and the 
waiting-time pdf is exponential, the fluid limit (in which only long-time, large-distance 
information is retained) yields diffusive equations. These choices intrinsically imply the 
existence of finite characteristic length and time scales associated to the underlying 
microscopic transport mechanism. 

However, there is a large number of physical, biological and social systems 0, BJ El E] 
in which the dominant transport mechanism lacks either a well-defined characteristic 
length scale or a well-defined characteristic timescale or both. Instead, the characteristic 
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transport scales diverge with the system size and/or lifespan. For instance, this happens 
whenever transport takes place via avalanches. Transport events as a result have a 
maximum size that is only limited by the system size L, and a characteristic size that 
diverges with some power of L. This situation has been encountered, to list just a 
few examples, while investigating the transport of both particles and energy out of 
magnetically-confined fusion plasmas [9], in the propagation of forest fires [TO] , 
earthquakes [TTJ, solar flares [12] and the transport of magnetic vortices in type-II 
superconductors [TB]. It is possible to construct CTRW models that lack characteristic 
scales by choosing p and ip from the family of stable Levy distributions [13]. The 
fluid limit is then expressed in terms of transport equations that contain fractional 
differential operators [15], which are essentially non-local and non-Markovian. There 
is evidence, both numerical and experimental, that models generated in this manner 
provide a reasonably effective description of transport for some of the aforementioned 
problems © Q21 Q2| . 

Although the formulation of CTRW models in finite systems is straightforward [HI 
[19], their fluid limit equations and associated propagators are much harder to obtain 
than in infinite domains. In this work we will study the formulation of CTRWs on a 
finite one-dimensional system with periodic boundary conditions (the unit circle) and 
work out its fluid limit in a rigorous way for an arbitrary stable Levy distribution. 

We envision that the formalism introduced in this paper might find application 
in systems in which particles or energy are amenable to trapping and/or avalanching 
processes along a periodic direction. A possible example might be the investigation of 
heat transport in single-wall carbon nanotube (SWNT) nanorings, which have gained 
recent attention due to their interesting transport properties [20], such as Aharonov- 
Bohm effects, magnetotransport or establishment of persistent currents. Some molecular 
dynamical simulations have found that linear SWNTs seem to exhibit heat conductivities 
divergent with the nanotube length, which suggests that a characteristic length scale 
might be lacking [21] . If these observations apply to nanorings as well, the periodic 
CTRW here presented might provide an adequate model for their effective description. 
Other possible applications might be to the modeling of some types of transport in 
turbulent plasmas confined in a (multi-)periodic system, such as those found in the 
(quasi-spherical) sun or in Earth-based (toroidal) magnetic confinement devices with 
interest for fusion energy production |22j. 

The paper is organized as follows. In Section [2] we derive the generalized master 
equation (GME) for a time-translational invariant, separable CTRW on the circle. In 
Section [3] we work out the fluid limit of such family of CTRWs introducing, in particular, 
the appropriate form of the Riemann-Liouville operators on the circle. The computation 
of the propagator of the fluid limit equations in the homogeneous case is presented in 
Section HI We discuss its asymptotic behaviour and give analytical solutions for a = 1,2 



in the Markovian case. Section [5] is devoted to the conclusions. Finally, |Appendix A 



and Appendix B give a brief survey on stable Levy distributions and fractional calculus, 
respectively. 
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2. CTRW on the circle. Generalized Master Equation 

We consider that the dynamics takes place on the unit circle, S l , parameterized by 
9 G [0,27r). Sometimes it is useful to view the unit circle as S 1 = R/27T, where the 
quotient is taken with respect to the equivalence relation x ~ x + 2ir for any x G R. 
Functions on S 1 can then be identified with the set of periodic functions on R with 
period 2tt. Explicitly, given / : S 1 — > R, its extension to a periodic function on R is 

/: R -> R 

(1) 

x ^ /(N), 

where [x] is the (unique) representative of the equivalence class of x which belongs to 
[0, 2n). Obviously, if x G [0, 2tt), f(x) = f(x). Finally, every function f E S 1 considered 
in this paper is assumed to have a Fourier series expansion, 

1 OO 
m=— oo 

with coefficients given by the formula: 

f m = f(9)e me d6. (3) 



o 



As mentioned above, we restrict ourselves to the case of separable CTRWs with 
time-translational invariance. This means that the CTRW is defined by two pdfs, p(A, 9) 
and iJ)(t, 9), called the step-size pdf and the waiting-time pdf, respectively. Here, A e R, 
r > and 9 G [0, 2ir). p gives the probability that a particle located at 9 jumps to 9 + A 
after having waited at 9 for a lapse of time r. Since A is unbounded, this may imply an 
arbitrary number of turns around the circle. Finally, p and ip verify the normalization 
conditions, 

/oo 
p(A,0)dA = l, Wg[0,2tt), (4) 
-oo 

POO 

/ ^(r,#)dr = 1, V0 G [0,2tt). (5) 
Jo 

We denote by n(9, t) the density of particles normalized to the total number of 
particles, i.e. L n(9,t)d9 = 1, V£. The GME for the time evolution of n(9,t) 
corresponding to this class of CTRWs can be derived as follows [231 HE]- First, note that 
the normalized density of particles is equal to the probability of finding one particle at 
9 at time t. This probability can be rewritten as: 

n{9,t)= [ Q{9,t') V {9,t-t')dt', (6) 
Jo 

where Q(9, t') gives the probability that the walker arrives at 9 at time t' < t, and r](9, r) 
gives the probability of staying at 9 for a lapse of time r. Clearly, one has: 

V(0,t) = l- r)dr. (7) 
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Then, we rewrite: 

oo 

Q(e,t) = Y t <y(e,t), (8) 

3=0 

where Qi(9,t) denotes the probability that a single walker arrives at 9 at time t by 
performing exactly j jumps. Obviously, for an arbitrary CTRW the function Qi($,t) 
satisfies the following recurrence relation: 

Q j (e,t)= J2 de ' dt'pie-e' + 2nm,e')^(t-t',e')Q j - 1 (e',t'), (9) 

m=-oo Jo Jo 



where the sum in front of the integral explicitly accounts for the fact that particles can 
arrive at 9 from 9' through jumps of length \9 — 9' + 27rm|, m G Z. Q°(9,t) = 5(9)5(t) 
provides the needed initial condition. Adding over all possible values of j, we find an 
equation for Q(9,t): 

00 i>2n rt 

Q(9,t)-5(8)5(t)= V / d9' dt' P {9-9' + 2TTm,9 , )iP{t-t',9')Q{9',t'). (10) 

We proceed now to carry out the previously described periodic extension to the 
real line. Trivially, p, -ip and Q also satisfy Eq. ( ITOl) . In R, we can make the change of 
variables 9 = 9' — 2nm and use the periodicity properties of p, tjj and Q, together with 
the additivity of the integral with respect to the interval of integration, and convert 
Eq. (HO]) into: 

/oo rt 
d6' dt'p(9-9',9')^(t-t',9')Q(9',t'). (11) 
-oo JO 

This equation is formally identical to Eq. (12) from Ref. [18], which appears in the 
derivation of the GME in an infinite system. The main difference is that Q(9,t), whilst 
positive, is no longer a probability density in K, since as a periodic function it is no longer 
integrable. This does not invalidate the procedure, because the convolution theorem still 
applies. Consequently, we can take the final result: Eq. (24) of Ref. [18]. Namely, 

rt /*oo i*t 

d t n(9,t)= / dt' d9'p(9-9',9')<i)(9',t-t r )n(9',t') - <P(9,t - t')n(9,t')dt' , (12) 



(13) 



o 



where is also a periodic function defined by: 



t) = c- 1 



C[fj(9,t)} 



and £.[•] stands for the Laplace transform with respect to t. 

Again, h(9, t) is not a probability density function in R even when positive because 
it is periodic and thus non-integrable. We need to restrict the integrals to [0, 2n) to 
recover the probabilistic interpretation. To do so, we simply reverse the manipulations 
which led us from (HOD to ([II]), recasting the GME flT2D as 

rt p27T rt 

d t n(6,t)= dt' d9'p(9-9',9')(j)(9',t-t f )n(9\t')- 0(9, t - t')n{0 ,t')dt' , (14) 
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with 



p(0,9') = p(0 + 2Tvm,9') 



(15) 



From this expression it is obvious that p is positive-definite and periodic in its first 
argument, p(0, 9') = P(2tt,9'), for any 9'. Furthermore, the normalization of p in K 
implies the normalization of p in S 1 , 



2- 



p(6,9')d9 = 1, We [0,2tt), 



(16) 



so that p is a probability distribution on the circle. Hence, the GME ( fl4l) is written in 
terms of functions on S l and makes manifest that the dynamics of n(9, t) is well-defined 
on the circle. It is worth pointing out that ( ITBT) expresses the fact that particles can 
wind around the circle an arbitrary number of times, and hence is a consequence of the 
non-trivial topology of the manifold on which the CTRW takes place. Actually, p shows 
up as the effective step-size pdf on the circle and is related to the original step-size pdf, 
p, through (fT5l) . This relation is in fact the ballooning transform that was introduced in 
the plasma physics literature in the late 70s to deal with PDEs formulated on toroidal 
geometries [23]. Figs. 1(a) and 1(b) compare p and p when p is a Gaussian and a 



Cauchy distribution (see Appendix A ), respectively. As expected, the effect of topology 
is much more relevant for Levy distributions with algebraic tails than for the Gaussian. 





(a) Gaussian 



(b) Cauchy 



Figure 1. Comparison of p (dashed) and its associated p (solid), (a) p is a symmetric 
Gaussian distribution (i.e. a = 2). (b) p is a symmetric Cauchy distribution (i.e. 
a = 1). In both subfigures, p has a = 1 and is centered at 9 = ir. 



3. Fluid limit of the GME. Fractional differential operators on the circle 

Since we will be especially interested in the fluid limit of the CTRW when p and/or ip 
are Levy distributions (which only have a finite number of finite integer moments), we 
cannot carry out the calculation by the standard method of expanding the GME using 
moments of these distributions. The appropriate procedure involves going to the Fourier 
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(for space) and Laplace (for time) domains. The fluid limit, which implies retaining 
only the dynamical information pertaining to the longer times and larger distances, 
corresponds to taking the limits k — > 0, s — > 0, where k and s stand respectively for the 
Fourier and Laplace variables. 

Since p is periodic with period 2ir, we can expand it in a Fourier series with respect 
to its first argument: 

oo 

P(^) = ^- E Pm(0')e- im6 - (17) 

Z7T z — ' 

m=— oo 

Recalling the definition of p in terms of p, the coefficients of this expansion are: 
p m {9') = / p(9,9')e im6 d9 = / p(0 + 2nm',e')e ime d9 

oo r2-n{m'+l) °° p2n(m'+l) 

= J2 / p(6,8')e im§ e- l2nmm 'd6 = / p(9,9')e im§ d9 

m'=-oo ^ 2lTrn ' m'=-oo J 2nm ' 

f oo 

p{9,9')e imS d9=p{m;9'). (18) 

-oo 

where in the third equality we have performed the change of variables 9 = 9 + 27rm' 
and in the fourth equality we have used that e~ l27Tmm ' = 1. Hence, the m-th coefficient 
in the Fourier expansion of p with respect to its first argument is exactly the value of 
the Fourier transform of p with respect to its first argument at k — m. This relation is 
characteristic of the aforementioned ballooning transform |24j . 

Following Ref. [25], we proceed by inserting this expansion in ([141) . 
1 °° ft r^TT 

d t n(9,t) = — J2 J dt 'J MW,t-l!)Prtffl^^n(6',lf) 



m ' = — OO 

t 



t-t')n(9,t')dt'. (19) 
Now, we expand both sides in Fourier series with respect to 9 and equate coefficients, 

d t n m (t) =-L T ^ d9e^ m -^ e f dt> d9'<j>(9' ,t - t')p m/ (9')e im ' e 'n{9',t') 
^ n i Jo Jo Jo 



27 

m'=— oo 

"2-K 



- / d9e imd / <f>{9,t-t')n(9,t')dt'. (20) 
Jo Jo 

Since e< m - m '^d9 = 2ir5 mm , and p m {9') = p{m; 9') we obtain: 

rt /*2tt 

d t n m (t) = / dtf / d9'(j)(e',t-t'){p(m;9')-l)e imff 'n(e , ,t'). (21) 



Let A(fc; 0) be the characteristic exponent of p(A, 9), i.e. 6*) = exp A(fc; 0). The 
fluid limit is defined by taking p(k;9) ~ 1 + A(k;9). Therefore, Eq. (|2ip in the fluid 
limit is written as 

d t n m {t)= d9'A(m;9')e ime ' dt'cj)(9',t-t')n(9',t'). (22) 
Jo Jo 
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At this point we need to make choices for the step-size and waiting-time pdfs. 
A natural choice is to take stable Levy distributions, which satisfy the central limit 
theorem [T3]. Thus, following Ref. [26], we will assume first that the waiting-time pdf is 
a positive extremal Levy distribution (see Appendix A). Performing a (non-analytic, in 
general) expansion of the Laplace transform of ip around s = and keeping the lowest 
order terms: 

4>(s) « 7 6(0,1], (23) 

where A 1 = cos(7T7 /2) if 7 G (0, 1) and Ay = 1 if 7 = 1, which corresponds to an 
exponential waiting-time pdf. The convolution in the second integral of (1221) in the fluid 
limit yields: 



C 



dt'4>(9',t-t')n(9',t') 



A 7 r" 7 s 1 - 7 n(^ / ,s) 



(24) 



where we have made use of (1A.3|) . Then, Eq. (1221) in Laplace space reads: 



sn m {s) - n m (0) = A^'V" 7 / A(m; 9')n{9'; s)e m9 'd9'. 



2tt 



(25) 



Now, multiplying both sides by s 7 using (IB. 21) and (1B.3|) and Laplace inverting: 



-, — n m (t) = A~,t 7 
dt 7 v ; 7 



2 77 



A(m; 9 f )n(9',t)e im9 'd9'. (26) 

As for the spatial part, we will consider the case in which A(k; 9) is the characteristic 
exponent of a general Levy distribution (see Appendix A) with parameters dependent 
on 9: 

' -u a {9)\k\ a [1 - i/?(0)sign(fc) tan(22i)] + ifi(6)k a G (0, 1) U (1, 2], 
-a(9)\k\ +in(6)k 



A(fc;0) 



(27) 



a 



restricting ourselves to = for a — 1. 

Let us work out the case a^l. Eq. (126]) becomes 



dt 7 



n m (t) 



A, 



2tt 



a a (9')\m\ 



1 — iP(9')siga(m) tan ( 



/7ra 

It 



AmB' 



n(9',t)d9 f 



2 k 



im^(9'y md ' n(9\t)d9' ). 



(2« 



or equivalently, 



dt 7 



r 7 L 



( 



\m\ a (a a n) m -i tan y 



) |m| Q sign(m) {j3a a n) m + im (nn) m (t).(29) 



The inverse Fourier transform of Eq. (|29|) is the desired fluid limit of the CTRW on 
the circle. The analogue of (121?|) in an infinite system is Fourier inverted by means of the 
Riemann-Liouville operators defined in Appendix B. The following theorem (see [27] for 
related results) yields the suitable generalization of the Riemann-Liouville operators on 
the circle. 
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Theorem: Let f be a sufficiently well-behaved function on S 1 , a G (0, 1) U (1,2), 
and define the operators 

1 r 2lT ( 6 — 6'\ 



*v m = *r m - (ar) J [ WK (i + «, i + ^) d* 



where o-D" 27i-£)a flre ^ g ^^ emann _^^ ot(t ,^/ e operators defined in W. q is the 
smallest integer greater than a and ( is the Hurwitz zeta function, 

C{s,a) = J2 T ^ » Re(s)>l, m + a ^ 0. (31) 

T/jen, 

(o^/) m = (-im) a /m, ( 27r P a /) m = H7™. (32) 

Proof: 

First, notice that it is enough to show that V a e ime = (im) a e im6 and ^V a e im9 = 
{—im) a e imd . Let us give the proof for oD a . By the definition of the Hurwitz zeta 
function and recalling that T(z + 1) = zT(z), 



= C " D ° /W + r(^)5 5ki (2 I r(m+l) + i-9') 1+ "-' d9 '' ^ 

Performing the change of variables 6 = 6' — 2n{m + 1), and using the periodicity of / 
we have: 



oc 



v *l r fm - 

^ d6i J (2vr(m + 1) + # - 6') 1+a "<i 

±* r" W + *S m+1 »4S= Urn £ r IB S. (34, 

Recalling the definition fIB.lj) of the Riemann-Liouville operators we immediately obtain 
that 

P Q /(^) = ^D a f{6). (35) 
Since _ 00 D a e im,e = (im) a e ime , the result follows. An analogous calculation shows that 

Let us go back to Eq. (12"5|) . Making use of the Theorem and the identities 

\m\ = — — , m sgn m = ; rr^ 5 (36) 

1 1 2cos(7ra/2) ' 1 1 6 v ' 2zsin(7ra/2) ' v ; 
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for a G (0, 1) U (1,2), we can Fourier-invert Eq. (1291) in terms of the operators defined 
in ([H: 



d 7 

dpi K J 



2tt cos(7ra/2) 

T 7 



[ Z>? ((1 + P)a a n) + ^Pf ((1 - (5)a a n)} (9, t) 



(37) 



which is the fluid equation corresponding to the CTRW on the circle. The cases a = 1 
(with (3 = 0) and a = 2 can be obtained, respectively, as the limits a — > 1 + and a — > 2" 
of the above expressions. 



4. Propagator of the fluid limit equations in the homogeneous case 

In this section we will particularize the above results to homogeneous systems and will 
compute the propagator (also known as the fundamental solution) of the fluid limit 
equations of the CTRW on the circle. That is, the solution with initial condition 



tnm) 



n(0,O) = J2 6 ( 6 - 21 

m=—oo 

sometimes called the Dirac comb. Its Fourier expansion is 

1 oo 



-imd 



(38 



(39) 



i.e. n m (t)\ t=0 = 1, Vm e Z. 

If the system is homogeneous Eq. (126!) takes the simple form 

d 7 

— n m (t) = AyT 7 A(m)n m (t), 
which after Laplace transforming allows to give the solution of the propagator as 



n(0,t) 



2tt 



,7-1 



_s 7 - A 7 r~TA(m) 



(t)e 



-imd 



(40) 



(41) 



We can make further analytical progress by using the so-called 1-parameter Mittag- 
Leffler functions, defined by 



E p (z) := 



j=0 



r(pi + i) ; 



(42) 



for p > 0. These functions provide a suitable generalization of the exponential function 
(indeed, note that E\(z) = exp(z)) which is useful to us because of the property: 

' (t) = E p (±at 7 ) . (43) 



c- 1 



sP =Fa_ 

Thus, Eq. ( )4T|) can be recast as: 



n 



^ oo 



7 ) e~ im61 



(44) 
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In the Markovian case, 7 = 1, we can produce more explicit expressions. By setting 
7 = 1 in ()44p we obtain: 

n(9,t) = j- e A{m)t/T e- im9 . (45) 

m=— 00 

In particular, if A(m) is the characteristic exponent of a stable Levy distribution with 

13 = 0, 

00 

t) = — V e (- o H"+iHiA e -™e )aG ( 0i 2]. (46) 



m=— oo 



In the limit t — > 00, n(8,t) — ► 1 /27r, as required by the condition of conservation of 
probability 

Given a symmetric Levy distribution (/? = /j, = 0), it is well-known that the fluid 
limit propagator in an infinite, homogeneous system is invariant under rescaling of the 
so-called self- similarity variable, % := x a /t 1 , and the propagator (with initial condition 
localized at the origin) at x = decays as t~ x l a . Things work quite differently in a 
periodic system, in which the notion of self-similiarity does not make sense. For the 
time evolution of the propagator on the circle {(3 = \i — 0) at 9 — 0, n(0,t), there exist 
two distinct regimes: 

- Long-time asymptotics: If t » r/a a it is clear that 

n(0 I t)~^(l + 2e-* a */ T ). (47) 

Thus n(0,t) decays exponentially for any a. 

- Intermediate asymptotics: For t << r/a a , we can approximate the infinite sum 
which defines n(0,t) by an integral, 

n(0,t) = — V e - CTQ H Q */^± / e - Q H^ dn r i/« (48) 

Hence, at short times the a-dependent algebraic decay characteristic of the propagator 
in infinite systems is recovered. 

In the next subsections we work out the special cases a = 2 and a — 1, for which 
we can provide analytical solutions for the propagator. 

4-1- Diffusive transport: a = 2 

For a = 2 the step-size pdf becomes a Gaussian distribution and (146!) yields: 

^ 00 

n (6,t) = — ^ exp(-m 2 a 2 t/T + im(nt/r-9)) , (49) 

m=— 00 

which can be written in terms of the third Jacobi theta function [28] : 

00 

$3(2, g ) = 1 + 2 y g m " cos(2mz), \q\ < 1, zGC. (50) 



m=l 
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n(6,t) = ^(U--8),e-^) (51) 



Namely, 



2tt \2 \ t 

which is the solution of the heat equation in a periodic system. 
4-2. Cauchy distribution: a = 1 

Let us consider the case a — 1, /3 — 0, which corresponds to a Cauchy step-size pdf. In 
this situation Eq. (1461) becomes: 

oo 

?7,( " 



27T 

m=— oo 



oo oo 



t) — ^ e -\m\at/T+im(fit/T-e) ^ ^2) 

We can split the sum as 

ra (fl 7 t) = — ( j[] e (-^/r +J (M/r-e)) m + e (-*t/T-iQ*/T-0)) m \ ; (53) 

\m=0 m=l / 

so that both terms are geometric series with ratios exp(— at/r + i(fit/r — 8)) and 
exp(— crt/r — i(jit/r — 9)), respectively. Then, 

oo 1 

e (-<rt/T+i(jj,t/T-6))m = i 

^ l-exp(-at/T + i(iMt/T-6)) ' 

oo 1 

E J-at/r-i(jd/r-e))m _ i i (<-a\ 
~ 1 - exp(-*t/T - iQd/r - 0)) { ] 

and we can give a closed expression for the propagator: 

. 1 sinh(crt/r) 

n(6,t) = —, — . , V 1 / . -r. 55 

v ' J 2ix cosh(at/r) - cos{fit/r -6) K ' 

It is instructive to work out the approximate form of the propagator for small t and in 
the vicinity of 9 = directly from Eq. (1551) . obtaining: 

n ^ *) w " < Tl YT^TTl Wr ^ 

7r (at fry + (fit/r — 0) z 

This formally coincides with the propagator in an infinite system. In particular, 
n(0,t) oc 1/t for small t, in agreement with the general result (148]) . 



5. Conclusions 



Although there are numerous results concerning CTRWs in infinite systems, much work 
still remains to be done on their finite-size formulation. In this paper we have performed 
a detailed study of time translationally invariant, spatially inhomogeneous, separable 
CTRWs on the circle, in order to investigate how the topology can affect the form of 
the Generalized Master Equation and the corresponding fluid limit equations. 

We have found that, in order to formulate the CTRW on the circle, a non-trivial 
transformation of the step-size pdf must be carried out which is strongly reminiscent 
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of the ballooning transform familiar to plasma physicists. During the derivation of 
the fluid limit of these periodic CTRWs, it was shown that the standard Riemann- 
Liouville operators need to be reformulated on the circle. We have provided original 
expressions for them. In the homogeneous, Markovian case, we have been able to 
compute the propagators of the resulting fluid limit equations, even in closed form 
for some simple cases. Similar calculations may be performed in the non-Markovian 
case, in a straightforward manner, by using Mittag-Lemer functions. Finally, we have 
also hinted how these propagators may prove useful when determining the characteristic 
exponents in experimental situations, by focusing on their short-time behaviour. We 
feel that the formalism developed herein will find an application in physical systems 
in which both periodicity and anomalous transport phenomena are encountered. As 
mentioned in the Introduction, this could be the case of systems like SWNT nanorings 
and gravitationally or magnetically confined plasmas. 
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Appendix A. Levy skew alpha-stable distributions 

The family of Levy skew alpha-stable distributions (or simply stable distributions, or 
Levy distributions) is parameterized by four real numbers a G (0, 2], (3 G [—1, 1], cr > 0, 
and fi G KL Their characteristic function (i.e. their Fourier transform) is given by 1 1 4J : 



According to the Generalized Central Limit Theorem ([291 [2]), stable distributions 
are the only possible distributions with a domain of attraction. The index of stability, 
a, is related to the asymptotic behaviour of S(a, j3, a, fi)(x) at large x: 



for a G (0, 2). For a = 2, S(2, (3, a, fi) is a Gaussian distribution. 

The skewness parameter, f3, measures the asymmetry of the distribution. If (3 = 0, 
S(a, 0, a, fi) is symmetric with respect to /i, whereas if (3 = ±1 the distribution is said to 
be extremal. Extremal distributions when < a < 1 are especially interesting because 
they are one-sided: they are only defined for x > if [3 — 1 and for x < if /3 = — 1. 
Their Laplace transform is given by (// = 0): 



S(a, (3, a, ij){k) 



exp {-a a \k\ a [l - i/3sign(£;) tan(^f)] + ifik) a ^ 1, 
exp (— a\k\ [l + z/5^sign(A;)ln|A;|] + z/i/c) a = 1. 



(A.l) 




(A.2) 




(A.3) 
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Appendix B. Fractional differential operators on the real line 



The Riemann-Liouville fractional differential operators are defined as [30] , [E] : 

»D»f •= r /W dx' (B 1) 

x r(m-a)dx m 7 ;r (x'-x) a " m+1 ' 

where T is the Euler Gamma function and m G Z is the ceiling of a (i.e. the smallest 
integer greater than or equal to a). When a — > — oo and 6 — > oo the action of the 
Riemman-Liouville operators in Fourier space satisfies the important property: 

where stands for the Fourier transform. This property justifies the name "fractional 
differential operators" because they generalize the notion of differentiation to non-integer 
order. 

Finally, we give the definition of the Caputo fractional differential operator of order 

7 e(o,i) j3U: 

d^ (t) - fo^) y d^ (t } (T^r ' ( ) 

for which 



(a)=a 7 £[/](«)-* 7 -7(0), (B.3) 



where £[•] denotes the Laplace transform. 
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